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Abstract 

High transverse momentum (pt) single non-photonic electrons which have been measured in the 
RHIC experiments come dominantly from heavy meson decay. The ratio of their pT spectra in pp 
and A A collisions {Raa{pt)) reveals the energy loss of heavy quarks in the environment created by 
A A collisions. Using a fixed coupling constant and the Debye mass (m,/) ~ gT) as infrared regulator 
perturbative QCD (pQCD) calculations are not able to reproduce the data, neither the energy loss 
nor the azimuthal {V2) distribution. Employing a running coupling constant and replacing the 
Debye mass by a more realistic hard thermal loop (HTL) calculation we find a substantial increase 
of the collisional energy loss which brings the V2{pt) distribution as well as RAAiPx) to values close 
to the experimental ones without excluding a contribution from radiative energy loss. 

PACS numbers: 12.38Mh 
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I. INTRODUCTION 



The spectra of mesons and baryons which contain hght flavors (u,d,s) only and which 
have been produced in ultrarelativistic heavy ion coUisions at the RHIC accelerator show a 
remarkable degree of thermalization. Hydrodynamical calculations reproduce quantitatively 
many of their dynamical properties and their multiplicity is well described in statistical 
model calculations. Statistical equilibrium, however, means loss of memory and therefore 
they are of limited use for the study of the properties of the matter which is created in the 
early phase of the reaction. 

Heavy quarks, on the contrary, do not come to an equilibrium with the surrounding matter 
and may therefore play an important role in the search for the properties of this matter. 
Produced in hard collisions, their initial momentum distribution can be directly inferred from 
pp collisions. The deviation of the measured heavy meson px distribution in AA collisions 
(divided by N^, the number of binary initial collisions) from that measured in pp collisions, 
is usually quantified as Raa = daAA/{Nc dp^) / {dapp/ dp^) . Raa is a direct measure of the 
interaction of the heavy quarks with the environment which is created in AA collisions. 
The same is true for the azimuthal distribution, da/d(f> oc (1 + 2vi ■ cos(0) + 2t>2 ■ cos(20)), 
where the V2 parameter is referred to as "elliptic flow", because at production no azimuthal 
direction is preferred. The observed flnite V2 value is therefore either due to interactions 
with light quarks and gluons or due to coalescence at the end of the deconflned phase when 
the heavy quarks are reshuffled into heavy mesons. 

In the RHIC experiments heavy mesons have not yet directly been measured. Both, the 
STAR [1] and the PHENIX [2] collaboration, observe single non-photonic electrons only. 
They have been created in the semileptonic decay of heavy mesons. Thus experimentally 
one cannot separate between charm and bottom hadrons. pQCD calculations in Fixed 
Order + Next to Leading Logarithm (FONLL) predict a ratio of aii^/acc = 7xl0~^ with the 
consequence that above pt > cross ~ 4GeV electrons from bottom mesons dominate the 
spectrum [3]. The uncertainty of this value is, however, considerable. The little known form 
of the electron spectrum from heavy meson decay and the little known ratio of heavy quark 
mesons to heavy quark baryons [4] add to this uncertainty. 

In order to understand the single non-photonic electron spectra one has to meet two 
challenges: One has to understand the interaction of a heavy quark with the environment. 
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produced in heavy ion collisions, and one has to understand how this environment changes 
as a function of time. In the past, several theoretical approaches [5-14]. have been advanced 
to meet these challenges. Almost all of them assume that in the heavy ion reaction a quark- 
gluon plasma (QGP) is created and that the time evolution of the heavy quark distribution 
function, f{p,t), in the QGP can be described by a Fokker-Planck approach 

= ^^[Mmm + ^Mp)fim- (i) 

In this approach the interaction of a heavy quark with the QGP is expressed by a drag 
(Ai) =< {p — p')i > and by a diffusion {Bij = | < {p — p') i{p — p') j >) coefficient calculated 
from the microscopic 2 — > 2 processes by 

1 f d^k f £k' f £p' 
< X > = — ' ' ' 



2E J (27r)32A; J {2TT)^2k' J {2txY2E' 
n,{k) X {2ny6'^'\p+k-p'-k')j Y. l-^^l' ^ ■ (2) 

p{p') and E = Pq {E' = p'^) are momentum and energy of the heavy quark before (after) 
the collision and k(k') is that of the colliding light quark or gluon. di is 4 for qQ and 2 
for gQ. n{k) is the thermal distribution of the light quarks or gluons which is usually taken 
as of Boltzmann type. Aii is the matrix element for the channel i, calculated using pQCD 
Born matrix elements. Up to now the calculations are limited to elastic collisions (Qq and 
Qg). The matrix elements for these channels can be found in ref. [12, 15]. They contain 2 
parameters which have to be fixed: the coupling constant and the infrared (IR) regulator to 
render the cross section infrared finite. Up to now all calculations have used a fixed coupling 
constant, albeit different numerical values. As IR regulator usually a Debye mass rriD has 
been employed which is assumed to be proportional to the thermal gluon mass mo = PgT 
with (3 around 1. 

The Fokker-Planck approaches differ in the way in which the surrounding matter is taken 
into account. The Texas A&M group [6-8] uses an expanding fireball whereas the other 
groups [5, 9, 10] use hydrodynamical calculations, with different equations of state, however. 

Despite of different choices for as and m/j and of the different models for the expansion 
of the QGP it is a common result of all of these approaches that they underpredict by far 
the modification of the heavy quark distribution due to the QGP. One has to multiply the 
pQCD cross sections artificially by a K factor of the order of K ^ 10 (which depends on 
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the choice of as and of the IR regulator) to obtain agreement with experimentally observed 
values for Raa{pt) and for V2{pt) [5, 9, 10]. 

One possibility to reduce the value of K has been advanced by van Hees et al. [6, 7] who 
assumed that heavy D-mesons can be formed in the plasma and decay thereafter isotropically. 
One has, however, to await more precise lattice results to see whether such a nonperturbative 
process is indeed possible. 

It is the purpose of this article to improve these models in three directions: 1) we replace 
the Fokker-Planck equation by a Boltzmann equation because the momentum transfer is not 
well parameterized by the first and second moment only. 2) we introduce a physical running 
coupling constant, fixed by the analysis of e~^e~ annihilation and of the r decay, in the pQCD 
matrix elements. 3) we replace the ad hoc parametrization of the infrared regulator by one 
which yields the same energy loss as the HTL energy loss calculations [16, 17]. We will 
show that with these new ingredients pQCD calculations yield a larger stopping of heavy 
quarks in matter and bring the results of the calculation close to the experimental values of 
Raa{pt) and of V2{pt)- 

We do not address here the radiative energy loss whose importance is highly debated 
[5, 18, 19] because detailed microscopic calculations are not at hand yet. They may easily 
count for the factor of two which remains for Raa{pt) between the data and the calculation 
which includes collisional energy loss only. This will be the topic of an upcoming publication. 



II. INFRARED REGULATOR 



In order to calculate the drag and diffusion coefficients (eq. 2) using pQCD Born matrix 
elements [12, 15] the gluon propagator in the t-channel has to be IR regulated by a screening 
mass /i 

Frequently the IR regulator is taken as the thermal gluon mass [20] 



^ = ^ = f (l + l"/)4.asT^«^ (4) 



where rij {N^) are the number of flavors (colors) and mjj is the Debye mass. The infrared 
regulator is, however, not very well determined on flrst principles. Therefore, in the actual 
calculations [6-8, 12] /i^ was taken in between f^fT^ and with g^g = Anas ■ The IR 
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regulator is one of the main sources of uncertainty for the determination of the cross section 
(and hence for the drag and the diffusion coefficient) and it is therefore useful to improve 
its determination by physical arguments. 

For QED Braaten and Thoma [21] have shown that in a medium with finite temperature 
the Born approximation is not appropriate for low momentum transfer It has to be 
replaced by a hard thermal loop (HTL) approach to the gluon propagator. At high \t\ we 
can use the bare gluon propagator (left hand side of eq. 3). This approach we call HTL + 
hard calculation. To calculate differential cross sections using hard thermal loops is beyond 
present possibilities but Braaten and Thoma have shown that in QED the energy loss can 
be calculate analytically in the HTL + hard approach. Our strategy is now the following: 
We assume that the gluon propagator can be written in the form 

(5) 



t-Kml{T) ' 

and determine the value of k, by requiring that a pQCD Born calculation with this gluon 
propagator gives the same energy loss as the HTL + hard approach. 

We first deal with the QED case where the underlying hypothesis g^T"^ -C is more 
likely to be satisfied and focus our attention on the t-channel which is the only one suffering 
from IR singularities and therefore decisive for the choice of n. For the HTL + hard approach 
we follow ref. [21, 22] where the collision of a muon with an electron is calculated. Let us 
consider the energy loss 

dE. 1 f d^k r d^k' f d^p' 



dx 2Ev J {2'n)nk J {2Tif2k' J (27r)32E' 

nF{k){l - npik')) X (27iYs^'\p+k-p'-k')^ ^ \M,e-.,'e'f u; , (6) 

where v is the velocity of the heavy muon, u = E — E' is the energy transfer in the 
collision and d = 4 is the overall spin degeneracy. The total energy loss is the sum of two 
contributions: 

1. At small momentum transfer |t| = \{p — p'Y\ < \t*\, where \t*\ is an intermediate scale 
chosen between g^T^ and T^, the hard thermal loop regulates the infrared singularity 
and we obtain ([21, 22]) 



dE^ 



dx 



|*|<|**1 487r g^T 



l^Ti^- (7) 
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2. At large \t\ (|t|max > |^| > |^*|) no infrared regulator is necessary and we arrive at 
([21, 22]) 



dE,. 



dx 



\t\>\f\ 



9^ 
487r 



In 



SET 



1 



3 
4 



C(2) 
C(2), 



(8) 



Adding the HTL (eq. 7) and the hard (eq. 8) part, the intermediate scale t* disappears and 
we arrive at [23] 



dE„ 



dx 



V->1 



9^ 
487r 



HTL+hard 



In 



48ET 



g'T' ^ 4 C(2) 

We compare now this result with that obtained by introducing an infrared regulated gluon 
propagator eq. 3. In Born approximation we obtain the cross section: 



C'(2) 



(9) 



dap 
~dt 



9^ 



vr s 



M' 



2\2 



+ 



,2 + 2 



(10) 



We evaluate here the energy loss for the whole t-interval t G [tmim 0] and obtain (for details 
we refer to the appendix A) 



dE„ 



dx 



cff 



487r 



In 



SET 



^ _ 3 _ C'(2) 



(11) 



Comparing the pQCD Born (eq.ll) with the HTL 
has to be 



4 C(2) 

hard result (eq. 9), we find that fi'^ 



fx 



— — — mt, = — — 

2e ^ 



K 



1.0,2, 

2e 



(12) 



6e 2e ^ 2e 
in order to obtain the same energy loss in QED. 

Because QED and QCD have a very similar HTL-propagator structure the above ap- 
proach remains valid for QCD as well provided that ^ 1 and that /i^ is replaced by 
eq. 4. In the QCD case there is, however, the complication that we are, for temperatures 
achieved at RHIC, at the best at the borderline of the the range of validity of the HTL 
approach, m|, ^ T^. As a consequence, the HTL+hard model - commonly used by many 
authors - is in fact not independent on the intermediate scale t*. To demonstrate this 
problem we start out as in QED. For small |t| we obtain 



dEr 



dx 



Cpas 



|t|<|t* 



X 



:i 



x'^Y -It" 



dt{-t) 



PL 



X^)p2 



with V being the velocity of the heavy quark Q and the spectral functions 
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1 



— $5 



T^L{x) 



and Piit, x) = 53 

TT 



t + Ut{x) 



(13) 



(14) 



n^, and Ht are the self-energies evaluated in the HTL approximation: 
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(15) 



For large |t| we obtain (see eq. 6) 



Q 



dx 



\t\>\t* 



E 



f d^k 


f d^k' f 


1 {27iy2kJ 


{2'Kf2k' J 



d^p' 



^Q{\t\-\t*\) 



n. 



{k){lTn,{k')) X (2iTYs^''\p+k-p'-k')-Y.\M, 



(16) 



Here the matrix elements include qQ qQ as well as gQ gQ collisions. In contradistinc- 
tion to the QED case the sum of both terms depends explicitly on the intermediate scale t* 
in the region [m|),T^], as seen in fig. 1. 



dE 

— (GeV/fm) 
dx 

0.4r 



^^^^^^ y2 

HTL +"^^^>«,» 

^ semi-hard *^ 


ml 


B.T. 
(HTL+hard) 
station. 


"~" - 

T=0.25GeV 

p=20GeV/c 

a,=0.2 
mfl=0.45GeV_ , 


semi 
hard 
(v^>J))- 

1 1 


^^hard 

tail (v=o) 

* llTL 

1 1 



dE 

— (GeV/fm) 
dx 

0,41- 



0.3 



0.2 



0.1 



0.01 0.02 



0.05 0.1 0.2 
\t'\ (GeV^) 



0.5 



Q-.(2;rT) 
t-Kml(T) 




0.05 



T=0.25GeV 
p=20GeV/c 

a,=0.2 
mfl=0.45GeV 



0.15 0.20 0.25 0.30 



FIG. 1: Left: the total energy loss in the HTL + hard approach as well as the different components 
for a given choice of parameters as a function of the intermediate scale t* . The full lines are the 
sum of the HTL (blue dotted) and hard/semi-hard part (dashed purple for i/^ = 0, dashed-dotted 
red for z^^ ~ 0.16m|)). Right: The total energy loss evaluated with Born cross sections and with 
the propagator, eq. 5, as a function of k. Only t-channel contribution has been considered here. 



There we display the two parts of the energy loss (eq. 13, blue dotted and eq. 16, purple 
dashed) as well as the sum of both (purple full). Clearly, the total energy loss becomes 
stationary with respect to the intermediate scale only for a value of ~ 0.4 GeV^ S> 
T2(= 0.0625 GeV^) and hence in a region where the HTL approach is not valid anymore. 
Mathematically, this is due to the appearance of terms cx O(^) and (x 0{j^) which are 
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neither small nor do they compensate. Physically, we are in a regime where the interaction is 
screened over a distance of the same order as the mean distance between QGP constituents, 
so that a large part of the "hard collisions" will be affected by the medium polarization as 
well. Our prescription to cure this problem is to add an IR regulator to the hard part 
(as /i^ in eq. 3). We dubbed this approach therefore 'semi-hard'. The HTL part remains 
unchanged. The value of is chosen in that way that for a wide range of temperatures 
and heavy-quark momenta the sum of the HTL and semi-hard energy loss is independent 
of for \t*\ < i.e. in the range where the HTL approximation holds. The red bold 
line in figure 1 shows this independence of the total energy loss on t* when the hard part is 
replaced by the semi-hard (red dashed dotted) approach for p = 20 GeV, T = 0.25 GeV and 
z/^ ^ O.lGm'jj. We will adopt this value of u"^ for the further calculations. 

If we compare the t channel energy loss calculated in the HTL + semi-hard approach 
(shaded area in figure 1) with that obtained within our pQCD Born approach (eq. 5) we 
find a value of n around 0.15. This value is close to that obtained in QED (eq.l2). It is 
considerably lower than those used up to now in the pQCD cross section calculation. This 
is our first seminal result. 

III. RUNNING COUPLING CONSTANT 

The constant coupling constant as is the other quantity which limits the predictive power 
of the present calculations. In the published calculations as was taken in between 0.2 [21] 
and 0.6 [12] leading to a difference of a factor 9 for the drag and for the diffusion coefficient. 

As has been observed by Dokshitzer [24] there exists the possibility to define a running 
coupling which stays finite in the infrared by writing observables as a product of an universal 
effective time-like coupling and a process dependent integral. An alternative approach is to 
define an effective coupling constant, acsiQ^), from the analysis of physical observables. 
Two different experiments, e+e~ annihilation [25] as well as non-strange hadronic decays of 
r leptons [26], have been used to determine the infrared behavior of acfr(Q^) • The resulting 
coupling constants are infrared finite and very similar. These effective couplings are all-order 
resummations of perturbation theory and include all non-perturbative effects. We extend 
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the parametrization of the time-hke sector, Ref. [24], to the space hke sector, which leads to 

« ^ «eff(g') = forg^^O, (17) 
I ^ — TT atn(L+/7rj 

with /3o = 11 — I n/, ri/ = 3, and L± = ln(±g^/A^) . In the space-hke sector we replace the 
propagator 



t t - ' 

where /i^ is an IR regulator which we will specify below. The coupling constant ttcfr is 
displayed in fig. 2 for 2 and for 3 flavors. It has already been argued in [27] that a running 
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FIG. 2: dependence of the running coupling constant. 



coupling constant leads to the disappearance of the logarithmic E dependence of the energy 
loss at large energies: 



— oc as SvrT^T^ln^ 
dx m 



D 



— oc as[ii )T 



(19) 



with an IR regulator /i^ = [l,2]m|,, where the Debye mass m^) is determined self- 
consistently according to 



rhliT) = Y (1 + ^ 4vr a(-m^(T)) 



(20) 



However, this ambiguity of the coefficient leads to a non negligible uncertainty in the energy 
loss. 

In this work, we determine the optimal infrared regulator using the same strategy as 
for the non-running case: we calibrate the energy loss to the one obtained in a generalized 
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"HTL + semi-hard" approach this time with a running couphng constant. For this purpose, 
we assume that the (squared) Debye mass for a fixed couphng constant appearing in the 
hard thermal loop terms (eq. 15), m|,(T) = (1 + ^)g'^T'^, can be replaced by m?jj{T^t) = 
(1 + ^)47racE{t)T'^ . As illustrated on fig. 3 (left, full purple line), also here the total energy 
loss depends on the intermediate scale |t*| in the domain of validity of the HTL approach, if 
we employ the HTL+hard approach. Only if we replace the hard by a semi-hard propagator 

«eff(t) acs{t) . s 

t ^ t - Xml{T,t) ' ^ ^ 

we may obtain an energy loss which is independent on the intermediate scale t*. The optimal 
choice is A ~ 0.11 (see fig. 3, left, bold red line). 



dE 

— (GeV/fm) 
dx 



dE 
— ( 
dx 



-(GeV/fm) 



2- 



1.5 - 



HTL + 
semi-hard 



1 - 



T=0.25GeV 
I mc=0.76GeV 
p=20GeV/c 



istation. 




semi-hard 



1.5 



ttetf (t) 
t — KiriQ 



0.5 



0.02 



0.05 



0.1 



0.2 



0.5 



ffeEf (f) 



t-K{l+ -^)47r a,s(t) 



1^*1 (GeV^) 



T=0.25GeV 
mfl=0.76GeV 
p=20GeV/c 
. I 



0.05 



0.10 



0.15 



0.20 



0.25 



0.3 



FIG. 3: Left: Same quantities as in fig. 1 for the case of a running aeff- Right: The total energy 
loss in pQCD Born approximation for two different infrared regulators as a function of k. The 
shaded area corresponds to the energy loss calculated in the HTL + semi-hard approach (left). 



Using this prescription, the energy loss in the t-channel is found to be ~ 1.3 — 1.4 GeV/fm 
i.e. ~ 6 times larger than the energy loss found with the same parameters for the non-running 
coupling constant. For \t*\ < T^, the HTL contribution becomes negligible and the energy 
loss is given by the semi-hard part only (which is IR-convergent). Therefore, the natural 
IR regulator /i^ for our effective Born pQCD approach (eq. 18) is //^ = Km%{T,t), with 
K ^ A ~ 0.11, i.e. exactly the propagator of the rhs of eq. 21. 

However, the same energy loss can be obtained if one uses the simpler propagator of eq. 
18 taking /i^ = K,mj^(T) 

(22) 



acfr(t) 



t - Km%{T) 
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with K ^ 0.2 and the Debye mass defined self-consistently according to eq. 20. This is 
shown on the right hand side of fig. 3 and leads to our choice fJ'Qco — 0.2 m|,(T) for the 
propagator defined by eq. 18. We will show later that with these values the drag coefficient 
and hence the energy loss differs only slightly between these two models in the (T, p) range 
of interest for ultrarelativistic heavy ion collisions. We note in passing that a similar energy 
loss has been obtained by Wick et al. [29] in a simpler model for light quarks. 

IV. RESULTS 

In order to evaluate the consequences of our new approach we compare the results with 
those obtained for other choices of coupling constants and infrared regulators. They are 
summarized in table 1. From A ^ F the parameterizations become increasingly realistic. 
For the results presented below we include the s and u channels as well. They do not require 
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TABLE I: Coupling constants and infrared regulators used in our calculations 

any IR regulator and the coupling constant have been chosen as a — > acfr('5 — ^'^) and 
a — s> acs_{u — m?) because s = m? and u = correspond to the maximal "softness" in 
these channels. 

A. Cross sections 

The cross sections ^ for the different parameterizations of table 1 are displayed in fig. 
4, left for quarks and right for gluons. It is evident that both, a running coupling constant 
and a lower IR regulator, increase the cross section at small t whereas the increase at high 
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FIG. 4: Effective cross section da/dt for the different models (see table 1). 
t is rather moderate, but nevertheless visible in the gQ reactions, due to the u-channel. 



B. Individual collisions and transport coefficients 

For many interpretations it is interesting to see how the quarks loose their energy when 
traversing a plasma of a given temperature. For this purpose we study the differential 
probability Pi{w,p) that a heavy quark with a momentum p in the rest system of the heat 
bath looses the energy w by colliding with a plasma particle of type i: 



The condition H > 0, where 



H = {47r)*E' {s-{E + k)^)t^+({2Ek-s + m^'^ 



-Ak'^p^ + 2w{k{s + ml) - E{s - ml))) t - w'^{s - m^)^] , (24) 

with s = ml + 2Ek{l — cos6'(fc, p)), determines not only the limits t± in eq. (23), it also 
constrains the integral over k. 

The probability Pi{w,p) for c-quark with p = 10 GeV in a plasma of the temperature of 
T = 400 MeV is displayed in fig. 5. On the left (right) side we see the probability for cq (eg) 
collisions. Negative values of w mean that the heavy quark gains energy in the collision. 
Due to the u-channel contribution eg collisions are more effective to transfer a large amount 
of energy. The large majority of the collisions yield only a small energy transfer. To show 
which collisions are most important for the total energy loss of the c-quark we display in 
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FIG. 5: Differential probability P{w) that a c-quarks with an initial momentum of p = 10 GeV/c 
looses the energy u; in a collision with a plasma particle in a plasma at T = 400 MeV, on the 
left hand side for collisions with quarks, on the right hand side for collisions with gluons. For the 
different curves see table 1. 
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FIG. 6: Differential energy loss wPqiw) 



dEa 



(left) and its normalized integral (right), both 



9*^'^/ ~ dtdw 

evaluated for a heavy quark with an initial momentum of p = lOGeV/c colliding with a quark. 
For the different curves see table 1. 

fig. 6 (left) (the absolute value of) wPq{w,p) for cq collisions, (eg collisions would exhibit 
a similar behavior). This quantity is directly related to the differential energy loss: 



(25) 



dxdw 

Collisions with a small energy transfer become dominant when a running coupling constant 
is employed. Fig. 6, right shows 

dE, 



/w 
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(26) 
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and displays that collisions with an energy transfer of w < 1 GeV contribute 70% to the 
total energy transfer in our new approach whereas in the standard model (B) they contribute 
25% only. 

In order to make our calculation comparable with other Fokker-Planck calculations we 
present in fig. 7 the drag coefficient A as a function of the Q-quark momentum p (left 
for c-quarks and right for b-quarks). The calculation for the two fixed coupling constants 
as =0.3 and a5'(27rT) do not yield different drag coefficients as long as the IR regulator is 
the same. Therefore we do not pursue model A. If one changes the IR regulator from the 
standard value, m|,, to that reproducing the HTL energy loss (k = 0.15) one observes an 
increase by a factor of 2. A running as (acfr) with a standard IR regulator increases the 
drag coefficient for low momenta where the small-t exchanges are more important. If the 
low t collisions are enhanced by both, a running as and a small IR regulator, we see an 
increase of the drag coefficient by a factor of ~ 5. The drag changes not substantially if the 
IR regulator is calculated with a running coupling constant - model F - as compared with 
model E and we therefore discard model F from further calculations. If as remains fixed the 
drag coefficient remains moderate for all infrared regulators, as it does for a running as and 
the Debye mass as infrared regulator. B-quarks show a similar behavior but their drag 

A(GeV/fm) A(GeV/fm) 
3 1- 3 1- 




p(GeV/c) p(GeV/c) 

FIG. 7: Drag coefficient A (left for charm quarks, right for bottom quarks) as a function of the 
heavy quark momentum p. We display A for a temperature of T = 400 MeV and for different 
combinations of coupling constants and infrared regulators as defined in table 1. For the different 
curves see table 1. 

coefficient is - due to their higher mass - around 30-40% smaller than that of the c-quarks. 
For a given plasma-lifetime evolution, we thus expect a smaller energy loss of b quarks, but 
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it is far from being negligible, especially in the most realistic models E and F. 

The drag coefficient depends strongly on the temperature. In fig. 8 we display that of a 
c-quark with a momentum of 10 GeV/c. As expected in our model, a hot plasma is much 
more effective to quench a fast quark than a cold one. 

A(GeV/fm) 




0.1 0.2 0.3 0.4 0.5 

T(GeV) 



FIG. 8: Temperature dependence of the drag coefficient A for c-quark with momentum^ = 
10 GeV/c. For the different curves see table 1. The arrows show how the drag coefficients change 
if one replaces m/) by the IR regulator determined by the HTL -|- semi-hard approach. 

C. Nucleus-Nucleus Collisions 

After having discussed single Qq and Qg collisions we investigate now the consequences of 
our approach for heavy quark observables in ultrarelativistic heavy ion collisions. To study 
the time evolution of the heavy quark in a plasma, usually a Fokker-Planck equation has 
been used. This approach has several shortcomings: a) The drag and diffusion coefficients, 
calculated by eq.2, do fulfill the Einstein relation only in leading logarithmic order E/T [5]. 
This is not sufficient to assure the thermalization of the heavy quark [28] . Either one has to 
impose the Einstein relation or the asymptotic heavy quark distribution is a Tsallis function 
and not a Boltzmann distribution, b) Being a small scattering angle approximation (or, in 
other words, containing the leading order term of T / Eq only) the approach brakes down if 
the momenta of the q(g) and of the Q are of the same order, i.e. in the region where V2 
becomes large, c) Even for large energies Eq first and second moment only (eq.2) are not 
a good approximation to the energy loss. It can be seen in fig.4 (right) that hard transfers 
are not excluded in the gluonic channel, due to the QCD-equivalent of the Compton effect. 
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Therefore, for the calculation presented here, we use a Boltzmann equation approach as 
in ref.[ll] in a test particle version. In coordinate space the initial distribution of the heavy 
quarks is given by a Glauber calculation. For the momentum space distribution as well as 
for the relative contribution of charmed and bottom quarks we use the pQCD results of 
[3]. In the E866 experiment at Fermi Lab [31] it has been observed that in pA collisions 
J/psi mesons have a larger transverse momentum as compared to pp collisions. This effect, 
called Cronin effect, can be parameterized as an increase of < > by 6^ ~ (0.2 GeV)^ per 
collision of the incident nucleon with one of the target nucleons. For most of the calculations 
we then convolute the initial transverse-momentum distribution of the heavy quark [3] with 
a Gaussian of r.m.s ^ncon{f±) 6o. In this parametrization ncou is taken as the mean number 
of soft collisions which the incoming nucleons has suffered prior to the formation of the QQ 
pair at transverse position f±. Future studies of D/B meson production at RHIC may allow 
to improve on this choice. 

In our approach we then follow the trajectories of the individual heavy quarks in the 
expanding plasma, described by the hydrodynamical model of Kolb and Heinz [10, 30]. We 
parameterize the temperature T(r, t) and the velocity u^{r,t) field of this model and use 
this parametrization in a finite time step method to calculate the collision rate F ( eq. 2 
with X = 1) ioi Q + g Q + g and Q + q Q + Q reactions ([12, 15]) and for the 
different parameterizations of the cross section. For a given interval of the (Bjorken) time 
At, we then generate the number of collisions according to a Poisson distribution of average 
FAr and perform these collisions individually. When a collision takes place we determine 
the final momentum of the heavy quark by taking randomly a scattering angle with a 
distribution given by the cross section at a given temperature. In this method no small 
angle approximations are necessary and we arrive by definition at a thermal distribution if 
we place the Q-quark in infinite matter at a given temperature. 

As the time-point of the hadronization of the plasma is not well determined, we explore 
here two options: a) Hadronization of heavy quarks into D(B) mesons when the expanding 
system enters the mixed phase and b) at the end of the mixed phase. In the latter option 
more collisions are possible and we expect therefore a larger quenching of heavy quarks. 
Also for the hadronization we apply two approaches which give slightly different meson 
momentum distributions: a) either we apply exclusively the fragmentation mechanism as in 
p-p [3] or b) we apply the fragmentation mechanism for high momentum quarks only whereas 
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at low momentum heavy mesons are formed by coalescence. For this purpose we define the 
probability distribution g that a heavy meson of momentum P is formed by coalescence of 
a heavy quark with momentum pq with a light quark as 

9{PJq) =PI d'qniq,T)fiq-pQ)5{P-pQ - q) , (27) 

where T) is the thermal momentum distribution of the light quarks at the moment of 
hadronization and f is the probability density that the heavy quark with a momentum pq 
forms a heavy meson with a light quark of momentum q. In the calculation we evaluate g in 
the fiuid rest frame and take f as a boosted Gaussian. j3 is chosen such that g is normalized 
to unity for pq = 0. Finally the heavy meson undergoes a weak decay and creates the single 
electrons which are observed in the detector. 

The results for Raa in central Au+Au collisions are compared to the experimental data 
in fig. 9. From top to bottom we show the results for the approaches B-E of table 1. On 
the left hand side we present the results for an hadronization at the beginning of the mixed 
phase, on the right hand side that for an hadronization at the end of the mixed phase. We 
observe that the additional interactions in the mixed phase reduce the artificial K factor, 
shown in the figure, with which the pQCD cross section has to be multiplied to describe the 
data. For some of the curves we present the results for 2 different values of K, in others we 
show the influence of the different approaches for fragmentation, "frag" means that heavy 
mesons are exclusively created by fragmentation, "coal. + frag." means that they are rather 
produced by coalescence at low momentum. It is evident that the different hadronization 
scenarios have little influence on the K factor which is necessary to describe the data. For a 
constant coupling constant and the Debye mass as IR regulator (model B) one has to employ 
K-factors of the order of 10-12. A smaller IR regulator (model C) or a running coupling 
constant (model D) reduce this K factor to values of 5-10, still much too large in order to 
render the calculation understandable. Only the combination of both, model E, brings the 
K factor close to an acceptable value of 1-2, leaving nevertheless still room for radiative 
energy loss. 

A very similar observation can be made for the minimum bias calculations which are 
compared with the experimental data in flg. 10. On the left hand side we display the results 
for model B, on the right hand side for model E. For a flxed coupling constant and the Debye 
mass as the infrared regulator we need, as for central collisions, a K factor of around 12, 
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FIG. 9: Comparison of the experimental and theoretical results for central Au+Au collisions. We 
display Raa of single non-photonic e~ as a function of the heavy quark momentum pT- The purple 
line shows Raa for e~ from B-meson and the red line that of D-meson decay for the K values 
indicated in the figure. The blue line is the sum of both. On the left hand side we assumed 
hadronization at the beginning of the mixed phase on the right hand side at the end of the mixed 
phase. From top to bottom we display the resu^s for the parameterizations B-E (see table 1). 



whereas for the model E the K factor is reduced to 1.5-2. Thus for central and minimum 
bias calculations the same K-factors have to be employed, a minimal requirement for the 
validity of this reaction scenario. 

^AAlept ^AAlept 




2468 10 q2468 10 

PT(Ge\/c) pT(Ge\/c) 



FIG. 10: Comparison of the experimental and theoretical results for minimum bias Au+Au colli- 
sions. We display Raa of single non-photonic e~ as a function of the heavy quark momentum px- 
The red line shows the e~ for D-meson decay, the purple line those for B-meson decay and the blue 
thick line the sum of both. Hadronization is assumed to take place at the end of the mixed phase. 
On the left hand side we display the results of model B, on the right hand side that of model E 
(see table 1). The applied K factors are given in the figure, the Cronin effect is taken into account. 

The Cronin effect changes the Raa value only for momenta between 1 and 3 GeV, as 
can be seen in fig. 11. It is therefore without any importance for the understanding of the 
Raa values at large pt but brings Raa much closer to the data in the pt range where the 
V2 values are large. 

We come now to the discussion of V2- To our knowledge, the present theories based on 
pQCD have not succeeded to describe simultaneously the experimental Raa and V2 results. 
As shown in 12, left, for model B neither the Cronin effect nor an augmentation of the K 
factor beyond the value needed to describe Raa increases V2 considerably. What helps is a 
larger interaction time, i.e. a late freeze out. This is shown in fig. 12, right, where we compare 
the V2 values for a hadronization at the beginning and at the end of the mixed phase. Using 
a fixed coupling constant the K-factors remain, however, large. If one combines a running 
as with a HTL -|- semi-hard infrared regulator one can reproduce V2{pt) using a K-factor 
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FIG. 11: Influence of the Cronin eff'ect on the Raa of single non-photonic e~ as a function of the 
heavy quark momentum pT- The red hne shows the e~ from D-meson decay, the purple line those 
from B-meson decay and the blue thick line the sum of both. We assumed hadronization at the 
end of the mixed phase. On the left hand side we display the results without, on the right hand 
side with the Cronin effect, both for the model E (see table 1). 
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FIG. 12: Dependence of the V2 of single non-photonic e~ as a function of the heavy quark momen- 
tum px on the Cronin effect and on the K factor (left) as well as on the freeze out density and on 
the infrared regular (right). All calculations are done with asi^nT). 



slightly larger than 2 and assuming a late freeze out, as can be seen in fig. 13. 

One could imagine that azimuthal correlations of non-photonic — pairs created in 
the decay of the heavy mesons whose heavy quarks have been created together may carry 
information on the energy loss mechanism. Many collisions with small momentum transfer 
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FIG. 13: V2 of single non-photonic e as a function of the heavy quark momentum for different 
freeze out energies using a running coupling constant and a smah infrared regulator (model E). 



may better conserve the original back to back correlations than few collisions with a large 
energy transfer. As displayed in fig. 14, this is not the case. Model A and model E give 
about the same azimuthal correlation. This means, on the other hand, that correlations are 
a quite robust observable to test this reaction scenario and to confront it with other ideas 
like the AdS-CFT approach [32]. 
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FIG. 14: Azimuthal correlation of e"*" — e~ non-photonic pairs as a function of the relative angle 
for model E and for both, central(left) and minimum bias (right), collisions; Q and Q are assumed 
to be produced back to back and the non-photonic e"*" — e~ background from uncorrelated pairs 
has been subtracted. 
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V. CONCLUSION AND OUTLOOK 



In conclusion, we have found that it is possible to reduce the uncertainties inherent in 
present day calculations of the energy loss and of the V2{pt) distribution of heavy quarks 
traversing a quark gluon plasma by 

a) determining the infrared regulator by the requirement that it reproduces the energy loss 
calculated in the hard thermal loop + semi-hard approach 

b) using an effective infrared safe physical coupling constant which describes other data like 
the gluon radiation in e~^e~ annihilation and the non strange decay of r leptons. 

Results of calculations in which these new features are employed come close to the ex- 
perimental data for Raa{pt) as well as for V2{pt)- The K factor required to reproduce the 
data is in between 1.5 and 2. Up to now a simultaneous description of Raa and V2 has 
not been possible even with large K-factors. That the K-factor is above one may be due 
to radiative processes which are not included here but it may also be due to the lack of 
a detailed knowledge of the different physical processes involved. They include the initial 
distribution of charm and bottom quarks, their hadronization and the role of heavy baryons. 

This observation has importance far beyond the physics of heavy mesons. Because the 
same running coupling constant and the same infrared regulator appear also in the cross 
section for light quarks we expect a similar energy loss for light quarks. Pions show indeed 
a very similar RAAiPr) distribution but baryons do not. The reason for this is unknown but 
if one follows the idea that they are formed by coalescence their formation mechanism may 
be rather different compared to that of heavy mesons. This conjecture is supported by their 
large V2 values. In addition, the large collective radial flow counteracts to the individual 
energy loss. To clarify the hadronization mechanism of light hadrons one has probably to 
wait until jet like hadrons and those created by the plasma hadronization can be separated, 
either by measuring correlations or by extending the detection range in momentum space in 
future LHC experiments. 

The observed enhanced cross section may also be of importance for the understanding 
of the fast equilibration observed in entrance channel of ultrarelativistic heavy ion collision 
where we do not have a heat bath like here but a momentum distribution given by the 
structure functions. There the typical momentum is, however, not far from that of the heat 
bath particles. 
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VI. APPENDIX 
A. HTL+hard 

As tke large transfer t will bring tke parton to a final state k' for wkick npik') << 1 we 
neglect tke factor 1 — np{k') for tke final state particle. We start from 

\t\>\t'\ 
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dx 
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np{k) I dt{—t)dF-j- 

i min 
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X 32^^ 
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wkere M is tke mass of tke muon and np is tke Fermi-Dirac distribution for a massless 
fermion. As 
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dE„ 
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(30) 



wkere s = + 2Ek{l — cos6{p, k)) and wkere tke integral is performed in principle 
over a domain suck tkat |tmin| ~ s > \t*\. For E ^ M ^ 1^*1^? one can nevertkeless 
argue on a pkysical basis tkat tkere is enougk "kardness" in almost every collision in order 
to fulfill tkis condition and tke domain in wkick tkis is not tke case becomes negligible. 
We will tkerefore integrate over tke wkole k space as tke integral converges. Introducing 
u = 1 — cos 9{p, k) e [0, 2], tke angular integral leads to 

Jo \ \t*\ A) 

, / AEk + M^ M2 , AEk + W^ 
A-K I In — h In 



\t* 



AEk 



M2 



(31) 
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Substituting the variable k hj x = k/T, we obtain the expression 
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Because E is assumed to be ^ M'^/T and because the integral is dominated by intermediate 
values of X (x ~ 1), one can neglect the last term in the integrand and take M = in the 
first term and one arrives at 
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eq. 7 of ref. [21]. 



B. effective IR regulator 



The t-integration of eq. 10 yields 
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and we obtain for the energy loss 
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where s = + 2ETxu. We first notice that in X (with a = '''^"^2^ ^ ) is maximal 

and bounded at a = (s = M"^) and then decreases like /i^/s oc fi^/ET for larger values 
of s. It then brings a contribution oc fi'^/ET that is subdominant at large energies. In this 
regime, s is ^ M"^ for most of the {u, x) integration domain, so that the 3rd term of X can 
be replaced by its asymptotic 1/4 value and |tmin| in the logarithm can be replaced by s. 
Therefore, 
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and one realizes that — 

ax 



cff 



is nothing but the hard contribution eq. 30 with e/i^ 



and — >■ + k, M^. We thus can read off the result directly from eq. 8: 
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and obtain eq. 11. 
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